INCOMPLETE SELF-SIMILAR BLOW-UP IN A SEMILINEAR 
FOURTH-ORDER REACTION-DIFFUSION EQUATION 



V.A. GALAKTIONOV 

Abstract. Blow-up behaviour for the 4th-order semilinear reaction-diffusion equation 

(0.1) u t = ~u xxxx + \u\ p ~ 1 u in MxM + , p > 1, 

is studied. For the classic semilinear heat equation from combustion theory 

(0.2) Ut = u xx + u p («>0), 

various blow-up patterns were investigated since 1970s, while the case of higher-order 
diffusion was studied much less. Blow-up self-similar solutions of (|0. 1 [) , 

U-(x,t) = (T-ty^f(y), y= (T _, )1/4 , 

are shown to admit global extensions for t > T in an analogous similarity form: 

u+(x,t) = (t-T)-p^F(y), y= ¥ ^ w . 

The continuity at t = T is preserved in the sense that 

U-(x,T~) = u+{x,T+) = Ci\x\~^ for all x € R \ {0} (d = const. ^ 0). 

This is in a striking difference with blow-up for ([0.2)1 . which is known to be always 
complete in the sense that the minimal (proper) extension beyond blow-up is u(x, t) = 
+oo for t > T . Difficult 4th-order dynamical systems for / and F are studied by a 
combination of various analytic, formal, and careful numerical methods. Other non- 
similarity patterns for (jO.ip with non-generic complete blow-up are also discussed. 



I. Introduction: self-similar blow-up patterns of higher-order 

REACTION-DIFFUSION EQUATIONS AND MAIN RESULTS 

1.1. Fourth-order RDE and blow-up. This paper continues the study began in [5l,[T3| 
[T8l [25] of blow-up patterns for the fourth-order reaction- diffusion equation (the RDE-4) 

(1.1) Ut = —u xxxx + in Rx R + , where p > 1. 

For applications of such higher- diffusion models, see surveys and references in [5l [25] . In 
general, higher-order semilinear parabolic equations arise in many physical applications 
such as thin film theory, convection-explosion theory, lubrication theory, flame and wave 
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propagation (the Kuramoto-Sivashinsky equation and the extended Fisher-Kolmogorov 
equation), phase transition at critical Lifschitz points, bi-stable systems and applications 
to structural mechanics; the effect of fourth-order terms on self-focusing problems in non- 
linear optics are also well-known in applied and mathematical literature. For a systematic 
treatment of extended KPPF-equations; see in Peletier-Troy [36J. 

The RDE-4 (11.11) has the bi-harmonic diffusion operator — and is a higher-order 
counterpart of classic semilinear heat equation from combustion theory [41] 

(1.2) ut = u xx + \u\ p ~ l u in Kx E+, 

with many well-known properties including various aspects of blow-up behaviour; see a 
number of well-known monographs [U [381 [22j [34J [33j [2H EE E7J . Surveys in [5] and in a 
more recent paper [18] contain necessary information concerning relations between these 
models (11.11) and (II. 2p and description of similarity and other blow-up patterns for (II. ip . 
Note that another related fourth-order one-dimensional semilinear parabolic equation 

(1.3) u t = -u xxxx - [(2 - (u x ) 2 )u x ] x -au + qe su , 

where a, q and s are positive constants obtained from physical parameters, occurs in the 
Semenov-Rayleigh-Benard problem [30], where the equation is derived in studying the 
interaction between natural convection and the explosion of an exothermically-reacting 
fluid confined between two isothermal horizontal plates. This is an evolution equation 
for the temperature fluctuations in the presence of natural convection, wall losses and 
chemistry. It can be considered as a formal combination of the equation derived in [27J 
(see also [B]) for the Rayleigh-Benard problem and of the Semenov-like energy balance 
[3"9~t [T2"] showing that natural convection and the explosion mechanism may reinforce each 
other; see more details on physics and mathematics of blow-up in [25]. In a special limit, 
(11 . 3 p reduces to the generalized Frank- Kamenetskii equation 

(1-4) u t = -u xxxx + e u , 

which plays a role of a natural extension of the classic Frank-Kamenetskii equation 

(1.5) u t = Au + e u , 

derived in solid fuel theory in the 1930s, [11] . 

Thus, similar to the second-order equations pair (II. 2p . (II. 5p . where both have an equal 
physical significance, we choose the RDE-4 (11.11) as a leading model, though some of the 
results are naturally applied to (11.41) . Note that (II .ip can be considered as a non-mass- 
conservative counterpart of the well-known limit unstable Cahn-Hilliard equation from 
phase transition, 

(1.6) u t = -u xxxx - (\u\ p ~ l u) xx in Ex R + , 

which is known to admit various families of blow-up solutions; see [10] for a long list of 
references. It is key that the mass-conservation for L 1 -solutions of (11.61) naturally demands 
existence of an extension after blow-up, which was studied in [T7J. This is not that 
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straightforward for the current model (11.11) . Indeed, (11.11) is related to the famous mass- 
conserving divergence Kuramoto-Sivashinsky equation from flame propagation theory 

(1.7) Ut = —u xxxx — u xx + uu x in Ix R + , 

which always admits global solutions, so no blow-up for (11.71) exists. 

As usual in PDE theory, blow-up means that, in the Cauchy problem (the CP) for 
(II. ip (or (11.21) ) with bounded smooth initial function uo(x), the classic bounded solution 
u = u(x,t) exists in R x (0,T), while 

(1.8) sup xeR jv \u(x, t)\ — > +oo as t —> T~ , 

where T G R + = (0,+oo), depending on given data u , is called the blow-up time of 
the solution u(x,t). It is also convenient to use the quite popular nowadays auxiliary 
classification from Hamilton [29], where Type I blow-up means the solutions satisfying, 
for some constant C > (depending on u), 

Type I: (T — t)~\u(x,t)\ < C as t — > T~ , and, otherwise 

(1.9) _^ 

Type II: lim sup^j.- (T — t) p- 1 sup x \u(x, t) \ = +oo 

(Type II also called slow blow-up in [29]). In R D theory, blow-up with the first dimen- 
sional estimate in (11.91) was usually called of self-similar rate, while Type II was referred 
to as fast and non self-similar; see [23] and [38] . 



1.2. Main results: on incomplete blow-up and extended global semigroup. Un- 
like the previous papers, we are now more interested in a possible extension of blow-up 
solutions beyond blow-up, i.e., for t > T. For (ll.2p . this is not possible: blow-up of any 
its solution u > is known to be complete, i.e., the proper (minimal) extension is 

(1.10) u(x,t) = +oo in Rx(T,+oo). 

In fact, (11.101) is guaranteed by the Maximum Principle; see [22] and [TBI Ch. 6,7] for 
main concepts and results of extended semigroup theory and further references. 
Our main goal is to justify, that, rather surprisingly for us, for the RDE-4 (11.11) . 



1.11) self-similar blow-up is incomplete and admits self-similar extension for t > T. 



Actually, (11.111) well corresponds to Leray's scenario of self-similar blow-up and similarity 
extension beyond, proposed by Jean Leray in 1934 for the Navier-Stokes equations in R 3 ; 
see [311 P- 245] for the precise formulation and [17] for recent discussions. For equation 
(II. ip . (11.111) has an important corollary: 



1.12) the number of self-similar extensions for t > T is not more than countable 



(actually, it is finite, but a proof exists in the analytic case p = 3, 5 ... only). This con- 
firms a plausible existence of an extended semigroup of unique global proper ( "minimal" ) 
solutions of the Cauchy problem for (11.11) . 

We will also present some convincing facts that self-similar blow-up is a generic (struc- 
turally stable) one for (11.11) . so that (11.111) implies a certain possibility of existence of a 
kind of extended semigroup theory of unique global proper solutions of (11.11) defined for 
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all t > 0. For (11.21) . this is well known; see [23] as a guide. For (ll.ip . this question will be 
also discussed but essentially remains open. 

For ( 11. 2p . the property ( 11. lip is revealed for the RDE-2 in M. N only, i.e., for the equation 

(1.13) u t = Au + u p , where p > p Soh = g±§, N > 2; 

see [23]. However, even for (II . 131) . such an incomplete blow-up is not generic in the sense 
that, for almost all (a. a.) solutions, blow-up is expected to be complete. This is not fully 
justified rigorously, though the results imply this at least in the radial geometry. 

Thus, according to (II. lip , the RDE-4 (11. ip exhibits more flexibility in adapting blow- 
up solutions for t > T, though the analysis becomes essentially more difficult. To justify 
(11.111) in Section [2] (blow-up patterns for t < T) and [3] (global patterns for t > T), we 
will use a variety of methods including some analytical ones, but the final conclusions will 
eventually depend on careful numerical experiments. This is an unavoidable feature of 
the study of the 4th-order dynamical systems to be derived, which seems do not admit a 
fully mathematically rigorous investigation. As a key issue, we claim that 



1.14) incomplete self-similar blow-up ( II. lip for (11.11) has a pure dimensional nature. 



Note that earlier, by a similar reason (11.141) . existence of extended solutions was obtained 
[TT] for the Cahn-Hilliard equation (11. 6p for p = 3: 

(1.15) u t = -u xxxx - (u 3 ) xx in R x R + , 

for which existence of a countable family of positive blow-up similarity solutions was 
established in [10]. However, since (11.151) is divergent and hence conservative by preserving 
the total mass of /^-solutions, existence of solutions beyond blow-up is more natural than 
for the non- conservative and non-divergent model fll.ip . 

In Section HI we present other blow-up patterns for ( 11. ip . introduced first in [13], which 
do not admit global extensions beyond blow-up, so belong to the case of complete blow-up. 
Fortunately, such complete blow-up is most plausibly non-generic for (11. ip . 

We hope that the conclusion (11.141) and others will be helpful for understanding singular- 
ity formation and extension concepts for 2mth-order nonlinear evolution PDEs including 

(1.16) u t = -(-/\) rn u+\u\ p - l u (m>2), 

which are definitely short of new ideas concerning construction of extended semigroups of 
unique global solutions. Note that some previous results in [51 [13] were already oriented to 
arbitrary m = 2,3, ... , though the present case of the even m = 2 in (11.11) will provide us 
with some surprises, especially for centre manifold patterns. In fact, dealing with (II. lip , 
we show how to extend a blow-up solution beyond singularity towards hence guaranteeing 
existence of a unique (this is most desirable, but not easy to prove) continuation. 

2. Self-similar blow-up patterns 
2.1. Blow-up similarity solutions. For convenience, we reduce the blow-up time to 



T = 0, 
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simply meaning that, by shifting in time, the Cauchy problem for (11.11) is considered in, 
say, (—1, +00) x R with initial data posed at t — — 1. 

Indeed, similarity blow-up is the simplest and most natural one for scaling invariant 
equations such as (II. II) . where the behaviour as t — > 0~ is given by a self-similar solution: 

(2.1) u-(x,t) = (-t)-£if{y), y=(zfr74, 
where f(y) ^ is a solution of the ODE: 

(2.2) A^(f) = -f i4) ~ hf -^ L if+\fr 1 f = m K, /(±oo) = 0. 
We recall that, for (11.21) . such nontrivial self-similar Type I blow-up (12.11) . 

(2.3) u_(x,t) = (-t)-i^/(y) l y = j^m =► r-i^'-^iZ+i/rv^o, 

is nonexistent. This was first proved in Ad'jutov-Lepin in 1984 pQ (see [20] for first 
applications of the nonexistence to blow-up evolution). For (j!.13p . nonexistence in the 
subcritical Sobolev range p < ^| was proved in Giga-Kohn |28j in 1985. 

But this is not the case for the RDE-4 (II. ip . Note that (12. 2p is a difficult ordinary 
differential equation (ODE) with the non-coercive, non-monotone, and non-potential op- 
erators, so the problem is not variational in any weighted L 2 -spaces. 

In what follows, by fo{y) we will denote the first monotone symmetric blow-up profile. 
We also deal with the second symmetric profile fi{y), which seems to be unstable, or, at 
least, less stable than f . 

The ODE (I2.2p was studied in [5] by a number of analytic-branching and numerical 
methods. It was shown that (12. 2p admits at least two different blow-up profiles with an 
algebraic decay at infinity. We also refer to [25, § 3] for further centre manifold-type 
arguments supporting this multiplicity result in a similar 4th-order blow-up problem. 
Without going into detail of such a study, we present a few illustrations only (mostly 
taken from [18]) and will address the essential dependence of similarity profiles f(y) on 
p. In Figure dj we present those pairs of solutions of (12. 2p for p = | and p = 2. All the 
profiles are symmetric (even), so satisfy the symmetry condition 

(2.4) f(0) = f" (0) = 0. 

No non-symmetric blow-up was detected in numerical experiments (though there is no 
proof that such ones are nonexistent: recall that "moving plane" and Aleksandrov's Re- 
flection Principle methods do not apply to (11.11) without the MP). Figure [2] shows similar 
two blow-up profiles for p = 5. 

2.2. Dimension of the "good" asymptotic bundle at infinity. To explain the na- 
ture of difficulties in proving existence of solutions of (I2.2p . let us describe the admissible 
behaviour for y 3> 1. It is not difficult to show that there exists a 2D bundle of such "good" 
asymptotics (see details in [9., § 3.3]): as y — > +00, 

(2.5) f(y) = [C iy ~p^ + ...] + [C 22 /~i ( H) e~ a °y 4/3 + ...], where a = 3 • 2~f 
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u t =-u xxxx + ' U ' P 1u: s y mme ^ r ' c blow-up similarity profiles for p=1 .5 u r =-u xxxx + l u l P 1u: * w0 symmetric blow-up similarity profiles for p=2 
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y y 

(a) p = |, N = 1 (b) p = 2, N = 1 

FIGURE 1. Two self-similar blow-up solutions of (j2.2[) : p = 1.5 (a) and p = 2 (b). 



u t =_u xxxx + ' U ' P u: tw0 ^low-up similarity profiles for p=5 




y 

FIGURE 2. Two self-similar blow-up solutions of (|2.2p for p = 5. 

and Ci and C*2 are arbitrary parameters. The same symmetric bundle exists as y— oo. The 
first square bracket with algebraic decay is connected with the dominated linear terms: 

(2.6) _ j/' y __!_/ + ... = /(y)=C7 iy -^ + .... 

In the second square bracket in (12. 5p . we see a typical WKBJ-type two-scale asymptotics 
in ODE theory. The constant a® in (12. 5p is obtained by substituting into the principal 
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part of (I2.2p the pure exponential term: 

(2.7) -/W- Jy/' + ... = 0, with /(y)~e^ 4/3 =► a 3 = -J (f) 3 , 

whence the unique (real) root — ao with a negative real part. We do not go into details 
of the asymptotic expansions here, since such calculus are well known in deriving optimal 
exponential estimates of the fundamental solutions of higher-order parabolic equations, 
which were first obtained in Evgrafov-Postnikov (1970) and Tintarev (1982); see Barbatis 
[2J [3] for key references and updated results. 

Note also that (12. 5p reminds a typical centre manifold structure of the origin {/ = 0} 
at y = oo: the first term in (12. 5p is a node bundle with algebraic decay, while the second 
one corresponds to "non-analytic" exponential bundle around any of algebraic curves. 

Thus, a dimensionally well-posed shooting is characterized as follows: 

(2.8) 



% uimeiisiuiiany weii-puseu siiuuimg is cnaiacieiizeu as iuiiuws. 
Shooting: using 2 parameters in (12.51) to satisfy 2 conditions (12.41) . 



In view of the analytic dependence of solutions of (12. 2p on the parameters Cip in the 
bundle (jZgfl the problem (Q, (T23D cannot have more than a countable set of solu- 
tions. Actually, the numerics and branching-homotopy approaches [5] confirm that in 
wide parameter ranges of p > 1 and N > 1, there exist not more than two solutions: 

(2.9) f (y) with {C w ( P ),C 20 (p)}, and f x {y) with {C 11 (p),C 21 (p)}. 

More precisely, existence of both the first (generic blow-up) profile fo(y) and the second 
one fi{y) are obtained by a "//-bifurcation" approach, when the ODE (12. 2p is replaced by 

(2.10) -/« - pyf - JL / + \f\P~if = 0, 

with a parameter \x > 0. Linearization about the constant equilibrium and y-scaling yield: 

(2.11) f = (p-l)-** + Y, y=(±)*z =► B*Y + ±Y + O(Y 2 ) = 0, 

where B* = — D\ — -r zD z is the well-known Hermite-type operator with the discrete 
spectrum cr(B*) = {— |, k = 0,1,2,} and eigenfunctions being generalized Hermite poly- 
nomials (I4.6p . [8]; see more details in Section H~2l It then follows from classic bifurcation 
theory [ID] that bifurcations occur in (12. lip when — j- gets on the spectrum of B*, i.e., 

(2.12) — j- = — | ^> 3 bifurcation points /Jk = \ for k = 2, 4, 6, ... , 

where by natural symmetry reasons we take into account even fc's only. It turns out 
that the bifurcations at /x 2 = \ and \i± = | are responsible for existence of f and f\ 
respectively; see [3 § 5-7] for details. For convenience, let us note that the first //-branch 
of solutions of (I2.10p originated at // 2 = f (a subcritical bifurcation) is strictly monotone 
decreasing for \i G (0, |) giving the fo at [A — \. The next /i-branch originated at ^ — -g (a 
supercritical bifurcation) is not monotone and is increasing on some interval [i E (|, /■**), 
so there is another non-zero element on it at // = j, which is precisely the f\. In general, 



This is rather plausible via standard trends of ODE theory [7], but difficult to prove. For positive 
solutions, it is true, and is always straightforward for odd p — 3, 5, where the nonlincarity is analytic. 
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the existence and multiplicity study of solutions of (12.21) in [5J [THl [25] is difficult and 
tricky, so here and later on we will need to essentially rely on careful numerical evidence 
to check the actual matching of the flows. 

The eventual similarity blow-up patterns (12. ip are characterized by their final time 
profiles: passing to the limit t — > 0~ in (12. ip and using the expansion (12.51) yields 

(2.13) if d{p) ± 0, then u_(z,*) -»• d\A~^ as * -»• °~ 
uniformly on any compact subset of R \ {0}. 

2.3. d = 0: final time profile as a measure. If d — in ( 12. 51) . i.e., /(y) has 
exponential decay at infinity, then the limit in (12.131) is different and there appears a 
measure in the data u-(x, 0~): in the sense of distributions, 

(2.14) Ci(p) = 0: -^E-6(x), t->0~; £_ = /|/| £ ^ < oo. 
Taking into account the sign of the solution, we have a measure in the "mass" sense: 

(2.15) dip) = : - e_«J(x), t - CP; e_ = / ^ 0. 

It is difficult to prove analytically that (12.141) actually takes place at some p = p§ > 1, 
so numerical methods have been used [18] to support this idea. Namely, the following 
exponent, for which (I2.14p holds, was detected: 

(2.16) p*(l) = 1.40... , for which E_ = 45.4244 

In Figure El we show such profile foiy) in the case (12.161) . accompanied by the second one 
fiiy), for which C\ ^ 0. 

2.4. On existence of similarity profiles: classification of blow-up and oscillatory 
bundles. We now prove existence of at least a single blow-up profile fiy) satisfying (12. 2p . 
(12 .4p . We perform shooting from y = +oo by using the 2D bundle (12. 5p to y — 0, where 
the symmetry condition (12 ,4p are posed (or to y = — oo, where the same bundle (I2.5P with 
y i— > — y takes place). By / = fiy;d-> d)i we denote the corresponding solution defined 
on some maximal interval 

(2.17) y G (yo,+oo), where y = y id,C 2 ) > -oo. 

If yoid, d) = — oo, then the corresponding solution /(y; d, d) is global and can repre- 
sents a proper blow-up profile (but not often, see below). Otherwise: 

(2.18) y id, d) > -oo f{y, d, d) -»• oo as y^>y£. 
Note that "oscillatory blow-up" for the ODE close to y = y^ : 

(2.19) / W = 1/1^/(1 + o(l)), 

where limsup/(y) = +oo and liminf fiy) = — oo as y — > y£ , is nonexistent. The proof 
is easy and follows by multiplying (I2.19p by /' and integrating between two extremum 
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Blow-up profile f Q (y) for P=P 5 =1 .40..., E_=45.4244..., e_=41 .7381 ... 

12 | 1 1 1 




5 10 15 20 



y 

FIGURE 3. Two self-similar blow-up solutions of (|2.2I) for p = ps = 1.40..., when 
El holds for / Q . 

points (2/1,2/2), where the former one y\ is chosen to be sufficiently close to the blow-up 
value i/q, whence the contradiction: 

o<Kn 2 (i/i)~-^i/r 1 (i/i)<o. 

We first study this set of blow-up solutions. These results are well understood for such 
fourth-order ODEs; see |26j . so we omit some details. 

Proposition 2.1. The set of blow-up solutions ( 12.191) is four-dimensional. 

Proof. The first parameter is yo G K.. Other are obtained from the principal part of the 
equation (12.191) describing blow-up via (12 .2p as y — > y$. We apply a standard perturbation 
argument to (12.191) . Omitting the o(l)-term and assuming that / > 0, we find its explicit 
solution 

(2.20) f (y)=A (y-y )-^, A'' 1 = $(-^), $(m) = m(m - l)(m - 2)(m - 3). 

For convenience, the graph of $(m) is shown in Figure HI Note that it is symmetric 
relative to mo = |, at which $(m) has a local maximum: 

(2.21) «&(§) = i. 
By linearization, / = /o + K, we get Euler's ODE: 

(2.22) (y - y ) 4 F( 4 ) = pAl 1 = PH-^). 
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Funsction 4>(m)=m(m-1)(m-1)(m-3) 




-10 12 3 4 



FIGURE 4. The graph of function 3>(m) in (|2.20p . and towards solutions of (|2.22p . 

It follows that the general solution is composed from the polynomial ones with the fol- 
lowing characteristic equation: 

(2.23) Y(y) = (y-y )™ => $(m) = p$(-^). 

Since the multiplier p > 1 in the last term in f!2.23|) and m = — ^rj is a solution if this 
"px" is omitted, this algebraic equation for m admits a unique positive solution m + > 3, 
a negative one m_ < — -^j, which is not acceptable by (I2.20p . ad two complex roots m 12 

with Remi j2 = | > 0. Therefore, the general solution of (12 . 1 9[) about the blow-up one 
(12.201) . for any fixed yo, has a 3D stable manifold. □ 

Thus, according to Proposition l2.lt the blow-up behaviour with a fixed sign (I2.18P (i.e., 
non-oscillatory) is generic for the ODE (12.21) . However, this 4D blow-up bundle together 
with the 2D bundle of good solutions (I2.5P as y — ► ±oo are not enough to justify the 
shooting procedure. Indeed, by a straightforward dimensional estimate, an extra bundle 
at infinity is missing. 

To introduce this new oscillatory bundle, we begin with the simpler ODE (12.191) . without 
the o(l)-term, and present in Figure [5] the results of shooting of a "separatrix" that lies 
between orbits, which blow-up to ±oo. Obviously, this separatrix is a periodic solution of 
this equation with a potential operator. Such variational problems are known to admit 
periodic solutions of arbitrary period. 

Thus, Figure fixed a bounded oscillatory (periodic) solution as y — > +oo. When we 
return to the original equation (12.21) . which is not variational, we still are able to detect 
a more complicated oscillatory structures at y — oo. Namely, these are generated by the 
principal terms in 

(2.24) / (4) = _i f2/+ | / | P -i / + ... as 
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Shooting periodic solution of f' 4 '=|f| p 1 f for p=5 




FIGURE 5. A periodic solution of = I/]* -1 / separatrix: p = 5. 



Similar to Figure [SI in Figure [6], we present the result of shooting (from y = — oo, which 
is the same by symmetry) of such oscillatory solutions of ( 12. 2ft for p — 5. It is easy to see 
that such oscillatory solutions have increasing amplitude of their oscillations as y — > oo, 
which, as above, is proved by multiplying (12.241) by /' and integrating over any interval 
2/1) ife) between two extrema. Figure [7] shows shooting of similar oscillatory structures at 
infinity for p — 7 (a) and p — 2 (b). It is not very difficult to prove that the set of such 
oscillatory orbits at infinity is ID and this well corresponds to the periodic one in Figure 
[5] depending on the single parameter being its arbitrary period. 

By C 2 (Ci) in Figure we denote the values of the second parameters C 2 such that, 
for a fixed C\ G 1R, the solutions f(y; C±, C 2 ) blow up to ±00 respectively. These values 
are necessary for shooting the symmetry conditions (12 .4p . 

Thus, overall, using two parameters Ci )2 in the bundle (I2.5P for y ^> 1 leads to a 
well-posed problem of a 2D-2D shooting: 



(2.25) 



find C\ 2 such that: 



yo(C 1 ,C 2 ) = -00, and 
no oscillatory behaviour as y 



-00. 



Concerning the actual proof of existence via shooting of at least a single blow-up pat- 
terns fo(y), by construction and oscillatory property of the equation (12.21) . we first claim 
that in view of continuity relative to the parameters, 

(2.26) for any d > 0, there exists C 2 *(Ci) e (C 2 (d), C 2 + (Ci)) such that f"(0) = 0. 

We next change C\ to prove that at this C 2 (Ci) the derivative f'(0) also changes sign. 
Indeed, one can see that 

(2.27) f'(0;C u C* 2 (C 1 ))>0 for C x < 1 and /'(0; C h C 2 *(Ci)) < for G x > 1. 
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Shooting oscillatory solution: p=5, N=1 , C,=1 , C 2 = 2.8956362303297 

4r 
3 - 
2 - 
1 - 

S <>•= 



-3 - 
-4 - 

-5 5 10 15 

y 

FIGURE 6. Shooting an oscillatory solution at infinity of (|2.2|) : p = 5. 



Schooting oscillatory solution: p=7, N=1, 0,=!, C 2 =2. 941 16376644039 Schooting oscillatory solution: p=2, N=1, C^l, C 2 =0.1131218908 




p = 2 (b). 

Actually, this means for such essentially different values of Ci, the solution f(y; C\, C|(Ci) 
has first oscillatory "humps" for y > and y < respectively. By continuity in Ci, (I2.27I) 
implies existence of a Cj" such that 

(2.28) /'(0;Cr,C 2 *(C 1 *)) = 0, 

which together with (I2.26|) induced the desired solution. Overall, the above geometric 
shooting well corresponds to that applied in the standard framework of classic ODE 
theory, so we do not treat this in greater detail. However, we must admit that proving 
analogously existence of the second solution fo(y) (detected earlier by not fully justified 
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arguments of homotopy and branching theory and confirmed numerically) is an open 
problem. A more difficult open problem is to show why the problem ( 12.251) does not 
admit non-symmetric (non-even) solutions f(y) (or does it?). 



3. Self-similar extensions beyond blow-up 

3.1. Global similarity solutions. Following [23] devoted to (11.131) and actually using 
Leray's blow-up scenario [31] , we suppose that the simplest way of extending of self-similar 
blow-up patterns for t > = T is using again similarity global patterns: 

(3.1) u+(x,t)=t-**F(y), y = ^r, 
where F is a solution of the ODE 

(3.2) A+(F) = -FW + \yF'+ {F^F = in R, F(±oo) = 0. 

In comparison with the blow-up one (12.21) . in (I2.2p two linear terms in the middle have 
changed their signs. We will show that this essentially changes the dimension of the 
asymptotic bundles and hence overall matching results. Note that exponentially decaying 
solutions of (13.21) were already studied in [19] (these are important for extensions of the 
blow-up of type (12.141) ). but now we need other types of solutions with algebraic decay. 
For convenience, we always impose the same symmetry conditions 

(3.3) F'(0) = F"'(0) = 0, 

and by F (y) we will denote those profiles, which can be considered for the role of an 
extension of the blow-up profiles fo. Then F (y) is not necessarily monotone for y > 
and even its uniqueness is rather questionable; see numerical analysis below. 

3.2. Dimensional analysis of matching for C\ ^ 0. We first apply to (13.21) a simple 
test such as (12. 7p to get the dimension of the exponentially decaying bundle: 

(3.4) F(y)~ e »y 4/3 =► a* = J (!) 3 , 

where the algebraic equation for a G C admits two roots with negative real parts: 

(3.5) a± = ao(— § ± i -v), where ao = 3-2~i. 

Therefore, instead of (|2.5|) . the bundle of such asymptotic orbits is 3D: as y — > +oo, 

+ {j / -!({S) e -?»'[ C - 2CO s(!m3^) + G 3 sin (aoS yi] +...}, 

where Gi and C3 are arbitrary parameters. We consider here the case G\ 0, and will 
treat the special one C\ = 0, with the data- measure (I2.14p . later on. 

The point is that the first constant C\ 7^ in ( 13. 6p is fixed by the blow-up limit (12. 13j) : 

(3.7) u + (x,0 + ) = d \x\~&. 
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Blow-up profile f Q (y) and its global extension F Q (y), p=2, N=1, 0^77.76.. 
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FIGURE 8. Blow-up profile fo(y) and the corresponding global one Fo(y) with 
the same C\ = 77.76... for p = 2. 



Overall, (12. 131) and (I3.7P will provide us with the necessary continuity of the unbounded 
self-similar solution at the blow-up time in the present case with C\ ^ 0: 



(3. 



) = u + (x,0 + ) for x 0, i.e., almost everywhere (a.e.) in 



Then {/o, -^o} is called a global extension similarity pair, or simply an extension pair. 
Thus, we again arrive at a dimensionally well-posed shooting: for a fixed values C\ ^ 0, 

(3.9) 



Shooting: using 2 parameters 6*2,3 m (I3.6I) to satisfy 2 conditions (I3.3I) . 



3.3. Numerical analysis of self-similar blow-up extension for C\ 7^ 0. However, 
the actual solvability of the problem lying behind (I3.8I) is difficult for a rigorous analytic 
study. The main point is that the solvability is very much Ci-dependent (and mostly is 
nonexistent for large Ci), so in many cases we will need again to rely on careful numerics. 

In Figure [H we show positive results of shooting the pair fo(y) and F Q (y) for p = 5. For 
comparison, we also put therein some other global profiles F(y) corresponding to other 
values of C\ in (I3.6I) from C\ = 1 up to C\ = 90, i.e., larger than the required 77.76... . In 
the next Figure [9l we show the enlarged asymptotic tails of all those profiles. 

Figure [TD] shows the extension pair {/o ; -^o} f° r V — 5, where C\ = 1.843... . In Figure 
[Til for the same case of p = 5, we present extra global F-profiles for smaller C\ = 0.2 and 
larger C\ = 2 values. We should note that convergence for the global F-problem is very 
slow and we succeeded in getting a few reliable numerics only. In all the cases, we have 
used the bvp4c solver of the MatLab with the enhanced accuracy with tolerances up to 

(3.10) Tols ~ 10" 10 . 
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x 10 -5Blow-up f (y) and global F (y), p=2, N=1, 0^77.76...; tails enlarged 




FIGURE 9. Enlarged tails of blow-up profile fo(y) and the corresponding global 
one Fq(jj) from Figure [8] with the same C\ = 77.76... for p = 2. 



f (y) and F (y) for N=1, p=5; 0=1.843... 




FIGURE 10. Blow-up profile fo(y) and the corresponding global one Fn(y) with 
the same C\ = 1.843... for p = 5. 



Note that not all the similarity blow-up profiles are assumed to have a global extension. 
For instance, Figure [T2l explains nonexistence of a global profile Fi(y) for the second blow- 
up one fi(y) from Figure [2] for p — 5. We present here the results of non-converging for 
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(a) profiles (b) enlarged tails 

FIGURE 11. The pair {f (y),F (y)} with C\ = 1.843... for p = 5 plus extra F-profiles. 



u t =-u xxxx + |u| p 1 u, p=5: blow-up ^(y) and non-converging global F(y) 



f(y) 



yy) 



F(y), non-converging 



0^4.446 



10 12 14 16 



FIGURE 12. Towards nonexistence of F\ for the second blow-up profile fi(y), 
with Ci = 4.446... for p = 5. 



the equation (13.21) with a sufficiently large 

d = 4.446.... 

Recall that, for G\ = 2 from Figure fTTT b). such a profile F (y) exists, but seems nonex- 
istent for larger values of C\. Possibly, this means that there exists a critical maximal 
value of C\ determining the optimal upper bound for existence of F-profiles. 
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3.4. On extension of a measure for C\ = 0. Consider the special case (I2.16p . where, 
as shown in (12. 141) . 

(3.11) p = Ps = 1-40... : lu-foO")!*^ = E-8(x), where E_ = 45.4244... . 

The corresponding blow-up profile fo(y) has been already shown in Figure [31 

If C\ — in (12.51) . one then needs to get the global extension profile F (y) also with 
C\ = in (13. 6p . This problem was studied in [19] by a bifurcation-branching approach. 
It was shown that there exists a countable sequence of critical exponents 

(3.12) Pl = l + £- v 1 = 0,1,2,..., 

such that in {p > pi} there exists a global p-branch of solutions Fi(y) of (13. 2\\ . that can 
be extended unboundently as p — > +oo. The behaviour of the branches near bifurcation 
points (I3.12p was proved to be: 

(3.13) F ol {y) = ll {p-p l ) 1 ¥ [^(y) + o(l)} as p -> p+; / = 0,1,2,..., 

where 7; > are some constants and V'j(y) are eigenfunctions (14. 9p of the operator (14. 7p . 

Before constructing the extension pair for the case (13.111) . note that this p$, within our 
accuracy, is close to the spectrum (13.121) . and, precisely, 

(3.14) p 9 = 1 + ± = 1.40. 

Note also that the critical exponents (13.121) are concentrated about p = 1 + for large /, 
and, in view of the asymptotics (13.131) . the p-branches are rather plain close to p = ps- 
Overall, this shows that there exist many global profiles Foi(y) satisfying (13.111) . 

Figure [TBI shows a typical global similarity profile F (y), though by the branching (I3.13p . 
there exists a countable set of such solutions of (13.21) with exponential decay. Therefore, it 
is difficult to identify which p-branch this profile F Q belongs to. The measure characteristic 
of this F (y) in Figure [TBI is: 

(3.15) E + = 13.2893.... 

The computations have been performed with the enhanced tolerances (I3.10p . since many 
profiles Fo are small enough due to (13.131) for p w p$. On the other hand, by the same 
branching, for I ^> 1 with pi more closer to 1 + , there exist Fq's with arbitrarily large E + \. 

Indeed, we see that the value (I3.15P for a global extension does not match the blow-up 
one in (I3.1ip : 

(3.16) P=Ps = 1-40... : E„>E + . 

Recall that, due to (I3.12p . there also exist global profiles F i such that E_ < E +i . However, 
all these discrepancies do not undermine a possibility of extension of the data (12.141) in 
a self-similar way given by (13.11) . An effective scenario of such a matching was obtained 
and confirmed in pjj] for the related limit unstable Cahn-Hilliard equation (I1.15p . which 
preserves the mass of the /^-solutions. Namely, it was shown JTOj § 4] that the rescaled 
blow-up solution can converge in the local topology (14.121) to a self-similar profile with a 
different mass (L 1 -norm for positive solutions). In practice, this means that close to blow- 
up time (as t — ► 0~), negative humps are created that eventually disappear at infinity, 
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y 

FIGURE 13. A typical global similarity profile i*b(y) f° r P = PS = 1-40... . 

and hence do not violate the matching. This "mass drift" is performed in the rescaled 
variables as t — > 0~ , so does not mean existence of the actual mass transferring mechanism 
in the x-variable that might be unavailable in the PDE. There is no a rigorous proof. 

We expect that a similar mechanism can be explored to neglect the mass "defect" , and, 
for t > 0, we may observe a global profile Foi(y), which minimizes the loss of the mass 
E- — E + i. The positive mass defect will then create two negative /positive humps that 
disappear at infinity in the rescaled variable y as t — > + . In this mass interpretation, we 
have to use the measure-like data obtained in the limit (|2.15p . Recall that the whole set 
of global profiles {Fqi} at p = p§ is countable, so that this minimization makes sense. It 
looks then like a certain "mass-discontinuity" is assumed at t = + , though this happens 
in the y = x/t 1 ^ 4 variable and hence the real discontinuity is not available in the x-one, 
which is prohibited for this parabolic flow. 

There are several open problems here, and the further analysis of the actual extension 
pair {/ , Fqi) based on enhanced PDE numerics is quite necessary and desirable to clarify 
this difficult transition blow-up singularity phenomenon. 

3.5. Towards extended semigroup theory. Thus, the above results, at least at a 
qualitative formal level, imply that self-similar blow-up for (11. ip (which is assumed to 
be generic and most structurally stable; see the next section) admits a proper extension 
beyond blow-up time. Moreover, we also can claim that, for C\ ^ 0: 

(3.17) there exists at most a finite number of extension pairs {/o, F }. 

For C\ = 0, i.e., for p = p$, due to the critical bifurcation exponents (13.121) . we expect a 
countable sequences of the extension pairs, but the discrepancy of their "masses" at t — 
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could reinforce extra evolution mechanism of choosing the right ones. This remains an 
open problem. 

Note that, in general, we cannot guarantee that a pair is unique. Nevertheless, in view 
of (13.171) . even in the case of finite multiplicity of the profiles Fq, there is a hope that 

(3.18) there exists a "minimal" extension pair {fo,Fo}, 

with F being "ordered" in a certain geometric-metric sense, e.g., F with the same C\ 
but with a most "less oscillatory" structure. Recall that fo is always monotone. 

This minimality property can play a key role in an attempt (seems rather naive) to 
constructing of extended semigroup theory for blow-up solutions of ( 11. lft by using standard 
ideas of parabolic regularizations via smooth solutions {u £ } (cf. [23] for (11.21) ). 

(3.19) u £ :u t = -u xxxx + ip £ (u), where <p e (u) = ij^pr^ -» \u\ p ~ x u as e -> + 

uniformly on compact subsets. Since <p £ (u) is globally Lipschitz continuous, the CP for 
(13.191) with the same data u has a unique global classical solution u £ (x, t). Then a proper 
(minimal) solution of (11.11) can be formally defined as a hypothetical limit 

(3.20) u{x, t) = lim e ^o u £ (x, t) in Rx(— l,oo), 

where the limit is understood in a pointwise sense, since u(x,t) is unbounded for t > 0. 
Moreover, no more restrictions on the topology of convergence in (I3.20p can be imposed, 
since, for general solutions, the blow-up set: 

(3.21) B[u ](t) = {x <E R : \u(x, t)\ = oo}, 

has an unknown structure for t > T (hopefully, of zero measure, which is not proved yet). 
Recall that even for much simpler nonlinear heat equations including f 1 1 . 1 3 [) . the limit 
extended semigroups are discontinuous in time in general; see [T6"t Ch. 7] for examples. 

Proving existence of the limit (13.201) (even along a subsequence) and checking how it is 
related to the minimal F -extension of the blow-up self-similar solution (12.11) are difficult 
open problems (possibly, non-solvable in general). Anyway, it is indeed surprising that 
there is still a slight hope that extended semigroup theory of blow-up solutions for (11.11) 
can be at least partially developed along the lines of that for the second-order parabolic 
PDEs such as (II. 2p [23l ITS] , where the Maximum Principle was always key. 

3.6. On generalizations to Mr. Finally, as a comment and an introduction to a future 
research, in Figure we show the extension pair {fo,Fo} for the equation in R , 

(3.22) u t = -A 2 u+\u\ p - 1 u, 

in the case p = 5 for N = 3 (a) and 6 (b). Note the clear difference in the geometry of 
the corresponding global F -profiles in (a) and (b). This is again a sign showing that the 
extension profiles Fo can be nonunique. 

The similarity solutions u± remain the same and the ODEs for radial patterns {fo, Fq} 
are easily obtained; see [18] for details on blow-up ones U-(x,t). The ODE problems for 
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blow-up Uy) and global F (y) for N=3, p=5; C 1= 2.328... Extension pair {f ,F Q } for p=5, N=6; ^=2.91 ... 




5 10 15 20 5 10 15 20 25 30 



y y 

(a) AT = 3, C x = 2.328... (b) N = 6, d = 2.91... 

FIGURE 14. The pair {fo(y),F (y)} for p = 5 and N = 3 (a) and N = 6 (b). 

/o and Fo in the radial geometry are as follows (now y stands for \y\ > 0): 

_ /( 4) _ M_A _ iN^M^l f» + (^Vz3) f _ I yf _ _I_ f + = Q) 

f(0) = /"'(0) = 0, f{y) = C iy -& + ... as y ^ oo; 

_ F (4) _ 2(7V-1) pm _ (N-l)(N-3) p„ , (jV-l)(jV-3) _p/ , 1 _p/ , _J_ _p , | F | P -1 F = g 

F'(0) = F"'(0) = 0, /(y) = C l2 r^ + ... as y -> oo. 

The asymptotic bundles (12.51) and (13. 6p remain analogous to the case N — 1. 

In Figure [151 we show the pairs {/o ; ^o} again for p = 5 in dimensions iV = 8 (a) and 
N = 9 (b). Finally, in Figure [16], the extension pair is shown for p = 5 and iV = 15. 
Analogously to N — 1, the convergence for the global profile F (y) is very slow, possibly 
in view of the multi-dimension of the bundle (13. 6p . However, for such profiles, we always 
observed a multiple convergence to such profiles from various initial data and/or varied 
accuracy and tolerances used in the bvp4c solver. 



4. NON SELF-SIMILAR "LINEARIZED" PATTERNS: COMPLETE BLOW-UP 

A countable set of such non-self-similar blow-up patterns for (II. ip were formally con- 
structed in [13]; see also [18] for extra details and necessary historical comments. Below, 
we briefly describe some necessary new features of these patterns to reveal the main rea- 
sons for their complete blow-up. For the RDE-4, there is no hope to get a full rigorous 
justification of existence of such non-self-similar blow-up scenarios, so we feel free to per- 
form a detailed formal construction. More advanced expansion and matching techniques 
in this direction can be found in [13] and also [H], where such a construction applied to 
non-singular absorption phenomena (regular flows with no blow-up), so a full mathemat- 
ical justification is available therein. 
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Extension pair {f D .F } for p=5, N=8; C 1= 3.22... Extension pair {f ,F Q } for p=5, N=9; C^S.38... 




5 10 15 20 25 30 35 5 10 15 20 25 30 



y y 

(a) N = 8, d = 2.22... (b) N = 9, d = 3.38... 

FIGURE 15. The pair {fo(y),F (y)} for p = 5 and iV = 8 (a) and N = 9 (b). 



Extension pair {f ,F Q } for p=5, N=15; C 1= 4.38... 




y 

FIGURE 16. The pair {/ , F } for p = 5, N = 15. 



4.1. Nonstationary rescaling. Dealing with non-self-similar blow-up, instead of (12. ip . 
we use the full similarity scaling: 

(4.1) u(x, t) = (-t)~r*v(y, t), y = — ^ , r = - ln(-t) -> +oo as t -> 0~. 
Then v(y,r) solves the following parabolic equation: 

(4.2) v T = A (v) = -v yyyy - \ yv y - ^ v + \ v \ p " 1 v in R x R + , 
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where A_ is the stationary operator in (12.21) . so that similarity profiles are simply sta- 
tionary solutions of (14.21) . 

4.2. Linearization and spectral properties. Performing the standard linearization 
about the constant equilibrium in the equation (14.21) yields a perturbed equation: 

(4.3) v = U + Y, U = (P ~ 1)^ Y T = (B* + I)Y + D(y), 

where D(F) = c^Y 2 + Co = |(p — 1)~, is a quadratic perturbation as Y — > and 

(4.4) W = -D$-\yD y in L% (R), p*(y) = e"^l 4/3 , ae(0,3-2-f), 
is the adjoint Hermite operator with some good spectral properties [8] : 

Lemma 4.1. B* : i/r»(R) — > L^*(R) a bounded linear operator with the spectrum 

(4.5) tr(B*) = {A fc = -f , fc = 0,l,2,...} ( = <r(B), B = + iy£> tf + ± /). 
Eigenf unctions ipl{y) are kth-order generalized Hermite polynomials: 

(4-6) V»:(y) = ^r[y* + EjK , ^»)V], fc = 0,1,2,..., 

one? i/ie subset {ipl} is complete in L^*(R). 

As usual, if {V'fc} is the adjoint (Riesz) basis of eigenfunctions of the adjoint operator 

(4.7) B = -Di + \yD v + \l in Lj(R), with p = i 

with the same spectrum (14.51) . the bi-orthonormality condition holds in L 2 (R): 

(4.8) (ij k ,ipi) = 6 M for any k, I. 
The generating formula for ip^. is as follows: 

(4.9) Mv) = ^ F(k) (y)i k = 0,1,2,..., 

where F(y) is the rescaled kernel of the fundamental solution b{x,t) of D t + D*, i.e., 

(4.10) b{x, t) = t~zF{y), y=^ k = -b xxxx , b{x, 0) = 8{x) =}► BF = 0, / F = 1. 

4.3. Why self-similar blow-up is expected to be generic. This can be connected 
with the spectral gaps in cr(B* + 1) = {1 — |} of the linearized operator in (14. 3p . Indeed, 
we observe here four positive unstable modes: 

(4.11) k = and 1, and k = 2 and 3. 

The first group is not taken into account since correspond to standard instabilities relative 
to perturbations of the time T = and the space x = of the blow-up point. 
The crucial is the actual unstable mode with 

(4.12) k = 2, where A 2 = \ and ij)*{y) = ±y 2 '. 

We claim that this unstable mode ensures convergence to a non-trivial similarity station- 
ary profile, which we used to denote by fo(y); see extra explanations in [25, § 3]. There 
is no still a proof of this intriguing fact on such a heteroclinic connection /* — > fo, but it 
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sounds rather reliably. Hence, the stable manifold of /o includes a part of the unstable 
manifold of the constant equilibrium f* — (p— l) -1 /^ -1 ), so that linearized patterns to 
be constructed must be less stable than /q, which then is expected to represent a most 
structurally stable (generic) blow-up pattern. 

Note that the third unstable mode in (14. lip with k = 3 has nothing to do with the 
second similarity profile j\. Indeed, this mode corresponds to odd (anti-symmetric) per- 
turbations, while fi(y) is an even function. The origin of existence of fx is more subtle 
and is explained in [5] by branching and homotopy theory. Moreover, according to [251 
§ 3.2], fi(y) is connected with the trivial equilibrium via a semi-stable centre manifold 
behaviour to be presented below. 

4.4. Inner expansion. Thus, in the Inner Region characterized by compact subsets in 
the similarity variable y, we assume a centre or a stable subspace behaviour as r — > +oo for 
the linearized operator B* + / (see [14] for rigorous details in a simpler related problem): 



centre: Y(y, r) = a(r)ipt(y) + w x (A 4 = -1, k = 4), 
stable: Y(y, r) = -Ce AfcT ^,(y) + w 1 (X k < -1, k > 4). 

As usual for the classic equation (II. 2p . we restrict to even values only, 

(4.14) k = 4,6,8,..., 

since existence of patterns with odd fc's, i.e., having a non-symmetric blow-up structures 
is rather suspicious. However, for higher-order diffusion, such patterns cannot be excluded 
entirely, though this is not our business here. 

For the centre subspace behaviour in (I4.13p . substituting the eigenfunctions expansion 
into equation (14 .3p yields the following coefficient: 

(4.15) a = a 2 7o + where 70 = c ((^) 2 , ipi) a(r) = + ... . 

It is crucial that (14.151) shows a clear semi-stable (saddle-node) structure of the equilibrium 
a = 0, which essentially depends on the sign of the coefficient 70; see below. 

Note that for the matching purposes, one needs to assume that (see details in [T3l 125]): 

(4.16) if ip* k {0) > 0, then 7o > and C > 0. 
Actually, for k = 4 it is calculated explicitly that 

(4.17) 70 = -^136^ < 0, 

so that the centre manifold patterns with the positive eigenfunction 

(4.18) v!(y) = ii(y 4 + 24) 



correspond to solutions that blow-up on finite interfaces, [251 § 3]. It seems reasonable 
that a full justification of such a behaviour can be done along the lines of classic invariant 
manifold theory (see e.g. [32]), though can be very difficult. 
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Overall, the whole variety of such asymptotics is characterized as follows 
(4.19) Y(y, t) = a(r)y k + ... , where a(r) = 



-J- + ... for k = 4, 

Ce AfcT + ... for k>A. 



4.5. Outer region: matching. The asymptotics (I4.13P is known [13] to admit matching 
with the Outer Region, being a Hamilton-Jacobi (H-J) one. More precisely, in the centre 
case with k = 4, according to 04.131) . (I4.15p . we introduce the outer variable and obtain 
from (14.21) the following perturbed H-J equation: 

(4.20) e = ^ =► v T = -\tvt-£ l v + \v\*- 1 v + i(iZvt-vxa). 

Passing to the limit as r — > +oo in such singularly perturbed PDEs is not easy at all even 
in the second-order case (see a number of various applications in [24J). Though currently 
not rigorously, we assume stabilization to the stationary solutions /(£) satisfying the 
unperturbed H-J equation: 

(4.21) _i£/'__I I / + |/|P-i/ = o in R. 

This is solved via characteristics, where we have to choose the solution satisfying ( 14. 19(1 : 

(4.22) /(e ) = (p _i)-^i_i. e 4 + ___ ag ^ Q 

Hence, integrating (14.211) with the condition (14.221) yields the unique H-J profile 

(4.23) /o(0 = /*(l + c4 4 )"^ T , where c=i-(p-l)A. 



Note that since 7o < according to (14.171) . the profile (14. 23(1 blows up as c^ A — > — 1; see 
analogous phenomena in [25]. For the 2mth-order PDE (I1.16P with odd m = 3,5, ... , we 
have 7o > and then (I4.23P can represent standard blow-up patterns, [13] . 

Analogously, for the stable behaviour for k > 4 in (14. 13(1 . we use the following change: 

(4.24) ^ = e ± iir T y v T = - v + {v^v - e~ V T v ^ m . 

Passage to the limit r — > +oo and matching with Inner Region are analogous. The H-J 
profile is as follows: 

(4.25) / c (0 = /* (l + c£ fc )~i^, where c = C (p - 1)5^ > 0. 

Overall, according to matching conditions ( 14. 19(1 and ( 14.241) . the whole set of such blow- 
up patterns consists of a countable set of different solutions for k = 6,8, ... (as we have 
mentioned, the case k = 4 works for m odd only, and then k = 2m). 

4.6. Final time profiles: the rout to complete blow-up. As a natural counterpart 
of the above asymptotic analysis, we introduce the final time profiles of solutions [13, § 4]: 

Centre subspace patterns: k = 4. This analysis is formal, since by (14.171) and (14.231) 
with c* < such patterns are not bounded. However, we keep the derivation formally 
assuming that case 70 > 0, which applies to equations (11.161) with odd m = 3,5,7,..., 
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[13]. Thus, marching of the outer region with the Outer Region II for \x\ > small, yields 
the following behaviour: 

(4.26) u(x,Q~) = C*\x\~^ \ \n\x\ |^(1 + o(l)), 

where > is a constant depending on the parameters p and 70 (assumed here to be 
positive) only and hence is independent of initial data uq. A more careful passing to the 
limit as t — > CT in (14. II) by using the limit profile (14.231) shows that 

(4.27) C* = 

Stable subspace patterns: k > 4. Such stable subspace patterns can be truly 
constructed for all m = 2, 3, ... . Matching with the Outer Region II yields, as x — > 0, 

(4.28) u(x,0-) = C4x\-&(l + o(l)) for A; = 6, 8, ... , 

where constants C* depend on initial function uq through the earlier constant C in (I4.19p . 

Comparing the self-similar profile (I2.13P with those in (I4.26P and (14.281) shows the actual 
origin of complete blow-up: the later ones contain stronger singularities as initial data at 
x = posed at t — 0. This is easier to see for the case k = 4 (actually, nonexistent by 
(14.171) ). where the profiles differ by the unbounded factor 

1 

(4.29) ~ I In |x| l"- 1 — > +00 as x — > 0. 

Surely, the same and in a stronger manner happens for k = 6, 8, ... . 

Bearing in mind our self-similar blow-up patterns (12.11) with the asymptotics (12. 5p . the 
unbounded factor in (14.291) actually means that, in a certain natural but formal sense, we 
have to look for a similarity pattern with 

(4.30) d = +00. 

Moreover, for existence of a regular extension for t > 0, this would demand existence of 
a global patters (13.11) . (13.61) also with the coefficient (14.301) . Of course, this is not possible 
that somehow reflects the nonexistence of a proper extension, meaning complete blow-up. 

The actual proof of nonexistence of a regular solution of (II. ip with data (I4.28P is 
not easy. This assumes proving nonexistence of a corresponding very singular solution 
for (II. ip . which is a well-known issue for second-order semilinear and quasilinear heat 
equations. However, such results are also known for higher-order nonlinear parabolic 
equations with various nonlinearities, where the analysis without the Maximum Principle 
becomes essentially more involved. We refer to functional methods in [2T], where necessary 
earlier references on the subject can be found. 

5. Final remark: on evolution completeness of blow-up patterns 

After introducing and discussing all the blow-up patterns, we finally are in a position 
to announce another important aspect of our study. Namely, we claim that the whole 
set of self-similar and linearized patterns introduced above are evolutionary complete in 
a natural sense; see [15]. This purely means that the self-similar blow-up patterns (12.21) 
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together with all the linearized ones in Section H] are expected to describe all possible 
types of blow-up that can occur for the PDE (11.11) . We connect this claim with the 
standard completeness of the eigenfunctions (14.61) in L 2 p ,. However, since by (14.111) we 
have thrown away two unstable patterns for k = 2 and k = 3 (the modes with k = 
and 1 do not occur in the evolution), one needs to replace these two ones. And this is 
done by including precisely two nonlinear self-similar patterns (12.11) with profiles fo and 
fx. Note that this may look irrelevant since the mode with k = 3 is spatially odd, while 
fi(y) was shown to be even. In fact, this is not that important, since we always mean a 
certain symmetrization of blow-up behaviour as t — > T~ (no proof is still available), so 
that including two "nonlinear eigenfunctions" /o,i(?/) i s expected to be enough to govern 
non-trivial self-similar blow-up. 

What is more crucial is that a standard centre manifold behaviour is prohibited by the 
sign in (I4.17p . Therefore, finally, according to evolution completeness arguments: 

two similarity profiles /o,i(l/) replace two modes, with k = 2 and 4. 

Here we observe a usual "dimension preservation" by including into the countable set of 
linearized patterns (minus two modes for k = 2,4) two nonlinear eigenfunctions fo,i{y)- 

Of course, another and obviously inevitable key aspect of the completeness speculation 
above is nonexistence of Type II blow-up in (I1.9P for (11.11) . Indeed, Type II blow-up 
patterns would destroy any completeness of the above functional set. Note that such 
Type II patterns were obtained for (13.221) for some large p in the supercritical Sobolev 
range 

P>Ps oh = ^ (iV>4); 
see [T5| § 5,6] for details. We have some reliable evidence that such patterns are nonex- 
istent in the subcritical range p < psob and, in particular, for N = 1, i.e., for (II. ip . The 
point is that, in the natural rescaled sense, such patterns on smaller compact subsets 
around the point {0, 0~} must be governed by a regular stationary solution W satisfying 

(5.1) - A 2 W +\W\ P ~ 1 W = 0, W(0) = 1. 

In the subcritical range, it is guaranteed that radial solutions of (15.11) are highly oscillatory 
[26] . which makes very unlikely to match such a behaviour with the outer region to get 
an acceptable blow-up pattern of (11.11) in the {x, t, w}-variables. We again refer to [HI 
§ 5], and must admit that a rigorous proof of such a nonexistence is absent. 
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